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Optimal Battery Charging for Damage Mitigation 
Technical Field 

This invention relates to model-based control methods in the charging and 
5 discharging of batteries and fuel cells to provide improved performance and longer life. 

Background of the Invention 

Nickel-hydrogen cells, as shown in Fig. 1, have been used in space applications for many 
years. The main advantages of nickel-hydrogen cells are high specific energy, long cycle life, long 

10 functional lifetime, and fairly robust chemistry. The long cycle life has been used to advantage in 
low earth orbit (LEO) applications where a satellite will orbit the earth every 90 minutes. Roughly 
35 minutes of this orbit is spent in the earth's shadow, during which time the satellite must use the 
batteries for power. The other 55 minutes of the cycle is spent in sunlight, during which time the 
satellite is powered by solar panels that must also recharge the batteries. In these applications, up to 

1 5 50000 charge-discharge cycles are common for nickel-hydrogen cells. 

Nickel-hydrogen cells are essentially a hybrid of standard cells and fuel cells. The positive 
electrode is nickel (NiOOH), and the negative electrode is hydrogen that diffuses into a porous 
platinum electrode. The platinum only serves to provide a reaction site for the hydrogen gas, and 
does not participate in the reaction. The entire cell consists of a large sealed pressure vessel that 

20 contains the hydrogen gas and the electrode stack. The electrode stack alternates pineapple-slice 
shaped disks of nickel electrode, separator (containing the electrolyte, KOH), platinum electrode, 
and a hydrogen diffusion screen. 

Given the very significant costs of powering satellites and space stations it is very desirable 
to extend the life of nickel-hydrogen cells as much as possible to reduce the cost of replacing the 

25 batteries in a space station or requiring the use of a new satellite when an old satellite could still 
function. Accordingly, there is a need in the art for a system to carefully control the re-charging of a 
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nickel-hydrogen cell and any other electrical storage devices in a manner to mitigate damage during 

charging and extend the life of the storage device. 



SUMMARY OF THE INVENTION 
5 It is thus an object of the present invention to provide optimal battery charging for 

damage mitigation. 

Another object of the present invention, which shall become apparent as the detailed 
description proceeds, is achieved by method for charging an electrical storage device so as to 
extend the life thereof, comprising: developing an essentialized cell model structure of the 
10 electrical storage device; determining model parameters for charge-discharge data of the 
structure; and determining charge-discharge behavior of the structure in a voltage-charge 
plane. 

Other aspects of the present invention are attained by an optimal re-charging 
controller for an electrical storage device, comprising a simulation processor receiving an 

15 input signal sent to the electrical storage device and a measured output signal generated by 
the electrical storage device, wherein the simulation processor models a plurality of states of 
the electrical storage device and generates an estimated output signal so that the controller 
mitigates damage to the electrical storage device. The simulation processor includes a 
feedback component which receives the measured output signal and the estimated output 

20 signal to generate a correction signal. And the processor includes an observer component 
which receives the correction signal and the input signal to generate the estimated output 
signal and an estimated internal state signal, wherein the estimated output signal converges a 
plurality of modeled dynamic states to corresponding states of the electrical storage device. 
Along with this there is a damage rate sensor that receives inputs signals from the electrical 

25 storage device and estimated states from the observer component to provide the control 
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* system with an instantaneous damage rate which is then used to modify the charging signal 

so that damage is reduced during charging. 

These and other objects of the present invention, as well as the advantages thereof 

over existing prior art forms, which will become apparent from the description to follow, are 

5 accomplished by the improvements hereinafter described and claimed. 



Brief Description of the Drawings 
For a complete understanding of the objects, techniques and structure of the 
invention, reference should be made to the following detailed description and accompanying 
10 drawings, wherein: 

Fig. 1 is a photographic representation of a typical Nilfe battery configuration; 
Fig. 2 is a schematic diagram of the dynamic components of an NiH 2 cell; 
Fig. 3 is a schematic diagram of typical voltage profiles inside a cell; 
Fig. 4 is a schematic diagram of an essentialized performance model of an NiBfe cell; 
15 Fig. 5 is a graphical representation of the charge-discharge cycle for the 500 th cycle 

of the cell according to the present invention; 

Fig. 6 is a graphical representation of the charge-discharge cycle for die 10000 th cycle 
of the cell according to the present invention; 

Fig. 7 is a graphical representation of the charge-discharge cycle for the 40000 th cycle 
20 of the cell according to the present invention; 

Fig. 8 is a schematic diagram of the solid phase relationships for an NiOOH 
electrode; 

Fig. 9 is a graphical representation of the energy distribution during the discharging 
portion of a typical charge-discharge cycle; 
25 Fig. 10 is a graphical representation of the energy distribution during the charging 

portion of a typical charge-discharge cycle; 
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J Fig. 11 is a graphical representation of the cycle life as a function of depth of 

discharge (DOD); 

Fig. 12 is an outline of a slow optimization approach; 

Figs. 13A-D are graphical representations of exemplary charging profiles from a slow 
5 optimizer, at constant current charging and optimal charging, wherein Fig. 13 A shows the 
charging profile as a function of charge, Fig. 13B shows the charging profile as a function of 
current, Fig. 13C shows the charging profile as a function of damage rate, and wherein Fig. 
13D shows the charging profile as a function of voltage; 

Fig. 14 is an outline of a fast optimization approach; 
10 Figs. 15A-F are graphical representations of exemplary charging profiles from 65 to 

100 percent fully charged, using a faster optimizer, wherein Fig. 15A shows the charging 
profile as a function of charge, Fig. 15B shows the charging profile as a fiinction of current, 
Fig. 15C shows the charging profile as a function of damage rate, and wherein Fig. 15D 
shows the charging profile as a function of normalized accumulated damage, Fig,15E shows 
15 the charging profile as a function of voltage, and Fig 15F compares charge and voltage; 

Figs. 16A-F are graphical representations of exemplary charging profiles from 60 to 
95 percent fully charged, using a faster optimizer, wherein Fig. 16A shows the charging 
profile as a function of charge, Fig. 16B shows the charging profile as a function of current, 
Fig. 16C shows the charging profile as a function of damage rate, and wherein Fig. 16D 
20 shows the charging profile as a function of normalized accumulated damage, Fig. 16E shows 
the charging profile as a fiinction of voltage, and Fig. 16F compares charge and voltage; 

Figs. 17A-F are graphical representations of exemplary charging profiles from life 
extending charging methods, wherein Fig. 17A shows the charging profile as a function of 
charge, Fig. 17B shows the charging profile as a function of current, Fig. 17C shows the 
25 charging profile as a function of damage rate, and wherein Fig. 17D shows the charging 
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profile as a function of normalized accumulated damage, Fig. 17E shows the charging profile 

as a function of voltage, and Fig. 17F compares charge and voltage; 

Fig. 18 is a schematic diagram of a real-time observer structure; and 

Fig. 19 is a schematic diagram of an advanced control system using an advanced real- 

5 time observer. 



BEST MODE FOR CARRYING OUT THE INVENTION 

In order to understand the concepts of the present invention, a detailed discussion 
will be provided for an exemplary electrical storage device - nickel-hydrogen cell - which is 
10 modeled and the specific parameters of the device determined for complete and full 
evaluation. It will be appreciated that the model selected for this particular type of electrical 
storage device are equally applicable to other electrical storage devices and as such the 
teachings of the present invention may find a wide range of applications in this field of 
technology. 



15 



20 



NiH 2 Cell Chemistry 

In normal charge-discharge operation of a nickel-hydrogen cell, the electrochemical 
reactions are as follows: 

discharge ^ 

positive nickel electrode: NiOOH+H 2 0 + e" Ni(OH) 2 + OH" (1) 



discharge ^ 



negative hydrogen electrode: £H 2 + OH~ H 2 0 + e~ . (2) 



charge 



During discharge at the positive electrode, the NiOOH (nickel oxy-hydroxide) on the electrode reacts 
with water in the electrolyte at the surface of the electrode as it accepts an electron from the external 
25 circuit through the nickel substrate. The products of this reaction are the Ni(OH>2 (nickel hydroxide) 
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•> on the electrode, and OH" (hydroxyl ion) which difiuses away from the from the electrode through 

the electrolyte in the separator. At the negative electrode, the hydrogen gas reacts with hydroxyl ions 

in the electrolyte to form water and release an electron in to the external circuit. This process is 

reversed during charging. The electrolyte, which is often 31% KOH solution in water, contains 

5 potassium ions, K + , and hydroxyl ions, OH", dissolved in water. Clearly the water and the hydroxyl 

ions are important in the reactions at the electrodes. It is usually assumed that the potassium ions do 

not take an active part in the cell reactions, however its exact role in the overall cell behavior is 

unclear. 

During charging, overcharge reactions can occur. These occur whenever the charging 
10 current drives the above reactions faster than the required mass transfer can happen. This can occur 
both when the nickel electrode is fully charged, and when the charging current is too large. Thus, in 
overcharge, the nickel and the hydrogen do not participate in the reactions, and the electrolyte is 
reacting electrochemically at both electrodes. The overcharge reactions are as follows: 

positive nickel electrode: 20H" overcWe > H 2 0+^0 2 + 2e" (3) 

15 negative hydrogen electrode: H 2 0+^0 2 + 2e" overcharge > 20H~ . (4) 

Consideration of these reactions would seem to indicate that overcharge is not a problem, but the 
formation of oxygen in a hydrogen environment is potentially very dangerous. 

In a collection of cells composing battery, cell reversal can sometimes occur when a single 
cell discharges completely before other cells. This cell then becomes a load on the other cells. If the 
20 cell is assumed to be nickel limited, the reversal reactions are as follows: 

positive nickel electrode: H 2 0 + e" f ™" af > OH~ + £H 2 (5) 



negative hydrogen electrode: ^H 2 + OH" reversal > H 2 0 + e" . (6) 

25 

It should be observed that the hydrogen gas that is generated at the positive electrode is consumed at 
the same rate at the negative electrode. Thus there is no net increase in hydrogen pressure or any net 
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change in electrolyte concentrations during reversal, and the cell can be operated continuously in this 



mode. 

Performance Model for NiH 2 Batteries 

5 Dynamic Modeling Background 

Control systems are often used to improve the performance of many systems and devices. 
Typically the quality of the control is related to the level of understanding of the system being 
controlled. This system knowledge is contained in the system dynamic model. Useful models for 
control design are dynamic, that is, represented by differential equations, and are of reasonably low 

10 order (requiring only a few differential equations) while still capturing the important dynamics of the 
specific system. A more complicated model, such as a computational fluid dynamic (CFD) model, 
will have thousands of differential equations (or dynamic states). This type of model is useful for a 
detailed understanding of battery phenomena and for cell design. Such a model, however, can 
contain states that are impractical to control, and can lead to a control system that is unnecessarily 

15 complex. Furthermore, a CFD model is often only represented as a computer code, and effectively 
no equations are available to the control system designer. On the other hand, a less complicated 
model, such as an empirical input-output relationship, may not contain any differential equations at 
all. This static model will usually not allow good control system design. Thus, to design a good 
control system, a low order dynamic model of the system to be controlled is most useful. This model 

20 should capture enough of the dynamic behavior of the system to be accurate over the frequency range 
of interest, while not being overly complicated. 

In the area of battery and fuel cell modeling, there would appear to be two common types of 
model, CFD-based and empirical. The highly accurate CFD-based models are based on first 
principles, and accurately model the diffusion of all of the active species through whatever media are 

25 present. Such models require a large number of electrochemical parameters, including things such as 
ion mobilities and reaction rate constants. Also required is specific geometrical construction data on 
a given cell. These models are highly accurate and provide good representations of cell dynamics as 
long as the cell configuration remains constant Unfortunately, as actual cells are charged and 

7 



WO 03/103071 PCTAJS03/17028 

discharged, their dynamic parameters do not remain constant, and these models typically do not 

provide good representations for electrochemical cells over many charge-discharge cycles. 

Additionally, the exact construction varies from cell to cell, and thus these models do not quite have 

the actual cell information. In terms of using these models for control, there are also a large number 

5 of states due to the representation of many species diffusing through the cell. Furthermore, the 

spatial lumping of such diffusive systems usually yields a dynamic system that is very stiff. Such 

systems are difficult to simulate, and rather than programming such a set of differential equations, an 

implicit simulation method is implemented directly on the diffusion equations in the computer code 

without the generation of any useful mathematical model on paper. Finally, these CFD-based models 

10 are expensive to build and change, while also being slow to build and use. 

At the other modeling extreme from the CFD-based models are the empirical models. These 
model typically have a specific algebraic structure, and the model parameters are identified from 
given data. These are usually simple models that may provide reasonable static operating point 
representations of a given cell as it charges or discharges. Other than having a built in integral to 

15 track the charge, these models are usually not dynamic. Consequently, they can be considered as 
static maps, such as that of a transistor or of a compressor. Although they can predict an operating 
point, they do not accurately represent dynamic changes. Thus, the empirical models are not 
adequate for control system design. 

A less common approach to battery and fuel cell modeling is what could be termed a hybrid 

20 model, or a performance model. Models of this type attempt to keep the important first principles 
dynamics, while not becoming so complicated as to not be readily useful. Their parameters are then 
tuned to match experimental data. This can be called essentializing the model. This model has a 
chosen structure based on an electrochemical understanding of cell behavior. Given this dynamic 
model structure, the parameters of the model are determined for a given cell from charge-discharge 

25 data. An important aspect of this model is the representation of the many diffusion processes internal 
to the cell by a fractional order system. Consequently the model is of low order (few dynamic states 
with some fractional) and can be used for many different cell chemistries. A model of this type is 
fast, understandable, and amenable to control design. Furthermore, it is possible to update the 
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parameters in the model as the cell changes through its charge-discharge cycles. Thus the hybrid 

modeling approach should provide a useful model structure for understanding cell changes through 

time, while also providing a useful dynamic model for control design. 



5 Overview of NiH 2 Cell Dynamics 

As in most electrochemical systems, the dynamics of NiH 2 cells are quite complex. This 
section presents an overview of the electrical and material behavior in a cell. Fig. 2 gives a generic 
structure of the major components of the cell dynamics in normal operation. Here the electrical and 
material contributions have been separated for clarity. At each electrode, a voltage is generated 

10 between the electrode and the electrolyte based on the concentrations of reactants. The reactants then 
diffuse into the reaction site or away from the reaction site. The hydrogen gas must diffuse into the 
platinum electrode from the ambient chamber. The nickel compounds may not diffuse, but it is not 
clear if the reaction front at the nickel electrode is diffuse or not. Also, the hydroxyl ions and the 
water must diffuse from one electrode to the other through the electrolyte. The electrolyte diffusion 

1 5 represents a major resistance to the flow of charge through the cell. An additional source of electrical 
resistance is the charge transfer resistance at the electrode surface, across which there is also a double 
later capacitance. A voltage profile inside a typical cell is shown in Fig. 3. There it is shown that 
voltage is produced between the electrode-electrolyte interfaces. It is shown later that these voltages 
are dependent upon species concentrations on either side and on the current passing through. 

20 Additionally there is always an ohmic voltage loss in the electrolyte due to current passing through. 
It should be noted that this ohmic loss always corresponds to heating. 



Model Essentialization 

To extract the essence of the dynamics from the complex model given in Fig. 3, several 
25 assumptions are made. The major assumption is that one electrode dominates the dynamics of the 
cell. In this case, a hydrogen rich environment, and the nickel electrode is dominant in the dynamic 
behavior. At this electrode, only the amount of charged material is assumed, c, is tracked, not the 
discharged material on the electrode. Also, all of the diffusive behavior into one fractional order 
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term, c di is lumped together which alleviates the need to model several diffusion channels in the cell. 

Finally, the voltage generated at the electrode will depend on the concentrations of both the stored 

material and the diffusing material at the electrode surface, and will be used as the terminal voltage. 

With these assumptions, the essential model is shown in Fig. 4. 

5 In the Fig. 4, the electrical behavior is shown in the right side of the figure, and the chemical 

species behavior is shown on the left side. On the electrical side, the voltage at the terminals, v, is 

due to the dependent voltage source that is a function of electrical current, z, concentration of charged 

material, c„ and concentration of diffusing material, c di at the electrode. On the chemical side, an 

electrical analogy of the material mass transfer is used. In the analogy, mass flow rate becomes 

10 current (or flux),y, and material concentration becomes voltage. Thus the two dependent current 
sources on the chemical side are mass flow rates, j s and j dy defined by the current on electrical side. 
The charged material is stored in a capacitor, C, that represents mass storage, with the addition of the 
self-discharge resistor, R sd , representing an orifice for the loss of charged material. The diffusing 
material must pass through a fiactional-order element, ks 9 , that is an agglomeration of all the various 

15 diffusion pathways in the actual cell. This fractional-order element separates the concentration of 
diffusing material at the electrode, c di from the concentration in the bulk of the electrolyte, c bu i k . 
The effective chemical reaction in the cell then becomes 

discharge ^ 

S + D + e " reaction products , (7) 

* charge 

20 where charge conservation is assumed in the essentialized model. Here S represents the reacting 
stored material and D represents the reacting diffusing material. The stored material in a NiH 2 cell is 
often associated with NiOOH, and the diffusing material with either H 2 0 or OH", but this is only for 
mental imagery and not inherent to the math model. The reaction products could also be tracked if it 
were necessary for the accuracy of the model, but they are not believed necessary. Based on this 

25 approximate cell reaction, during discharge, the stored material reacts with the diffusing material and 
an electron. Thus, when current is leaving the battery (electrons entering), both the stored and 
diffusing material are being deleted at the active electrode. Note that it could have been assumed the 
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diffusing material was building up outside the electrode during discharge, rather than going to zero, 

but that would have added complexity. And it will be appreciated that a thermal model is not 

developed directly for the essentialized model. It is assumed that the heating effects are related to 

the cell damage, and that the damage model developed later will contain the necessary information. 

5 Heating is also considered later where energy flow, or power, in the cell is discussed. 

Electrode Behavior: Faraday's Law 

For the discussion of the next two sections, consider the electrode process : 

. * +b> < > Cc + (8) 

10 

This states that chemical species A and B are converted to species C and D plus electrons, at rates 
related to the coefficients a,b,c,d,n (it is also assumed that these compounds are appropriately ionized 
so that charge is conserved). It is important to understand that this equation represents both a 
chemical process and an electrical process. In fact, the reaction rates can be completely determined 

15 using the electrical process by recognizing that the chemical conversion can only occur if electrons 
are either arriving or leaving. Thus, the chemical conversion rates are controlled by, or measured by, 
the electrical current passing through a given electrode. All the rates of chemical conversion are then 
known via the following observation. For every n electrons leaving the electrode, c molecules of C 
are produced, d molecules of D are produced, a molecules of A are consumed, and b molecules of B 

20 are consumed. Then recognizing that the rate of electron production is related to electrical current i, 
the following rate equations result: 

_ . g d <l e - = _ d dD _ _ c dC _ + a dA = + bdB (9) 
dt n dt n dt n dt n dt 9 

Where q e . is the charge of a single electron. It should be noted that concentrations of chemical 
25 species are often used. Concentration is defined as moles of the species per unit volume. 
Concentrations are also often represented by [A] or by c A , which would represent the concentration 
of species A. Using concentrations, the previous equation becomes: 
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- im nF^--d^--c^- + a^- + b*± (10) 
dt dt dt dt dt 



where F is the Faraday constant relating the moles of electrons in a coulomb, Z^=96484.6 CoulAnole. 

5 Electrode Behavior: Electrode Equation 

Once the concentration rates on either side of the electrode interface are known, it is 
important to obtain an expression of the voltage in terms of these concentrations. This can be done 
by considering conservation of charge through the interface. Again assuming the reaction of 
Equation (8), the required relationship is determined by considering the fluxes,y, of the species in the 
1 0 forward and reverse reactions at the electrode, 

i = Jf-jr- CI) 

Assuming that the species fluxes are proportional to concentrations, yields 
15 i = k f c a A c b B -k r c c c c d D . (12) 

Here, the rate constants, k f can be related to the electrical potential across the electrode-electrolyte 
interface using free energy considerations, 

20 . k f = k 0 e^ F,RT ^ ) ^ ) (13) 

k r = v^'^'"^ ( 14 ) 

where R is the universal gas constant, T is the temperature in K, a is an electrode energy barrier 
symmetry factor, v - v 0 is die voltage deviation from some nominal voltage v 0 , and A© is a rate factor. 
25 Inserting these gives the electrode equation, 
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' = *o {c' A c h B c <"* r ><»-»><~*> - c e c c d D e iF,Kni - aHv -^) . (15) 



This equation is important to both the electrical and chemical behavior of the electrode-electrolyte 
interface as it relates the voltage across the interface to the concentrations of active species on either 
5 side of the interface and to the electrical current passing through the interface. It is often called the 
cuirent-potential characteristic, and the approach used is often referred to as the Butler-Volmer 
approach. 

Thus, for an electrode-electrolyte interface satisfying a reaction given by Equation (8), the 
required equations relating electrical current through the interface, electric potential across the 
10 interface, and species concentrations on either side of the interface are Equation (10) and Equation 
(15). These equations can also be used in a multiport configuration to relate and separate the 
electrical variables from the chemical variables. Consequently, these two equations are fundamental 
to representing the behavior of any electrified interface. 

1 5 Linear-in-the-parameters Electrode Equation 

The electrode equation is highly nonlinear, particularly when solving for the voltage in terms 
of the concentrations and the current It is almost always solved using numerical methods, such as 
Newton -Raphson iteration. As an alternative to doing this, we propose a linear-in-the-parameters 
approximate solution to the electrode equation. Thus, the form we have chosen is: 

20 

v = fc, + k 2 ln(l + 1 i |)sgn(0 + k 3 \n(c d ) + k A ln(l - c s ) (16) 

where the variables i 9 c si c d correspond to those in our essential model and the k*s are parameters to 
be determined from data. What this equation assumes is that the terminal voltage v, is equal to the 
25 sum of four terms, and is effectively a compromise between the Tafel and the Nernst solutions of the 
electrode equation. The Nernst solution assumes that the current is so small as to be negligible. The 
electrode equation is then easily solved for the voltage in terms of the concentrations, 
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Alternatively, the Tafel solution assumes that the current is large in one direction or the other, which 
means that one of the two exponential terms is negligible. With one of the exponentials equal to 
5 zero, the electrode equation is easily solved for the terminal voltage in terms of the current, 

< I8 > 

In the linear-in-the-parameters approximate solution to the electrode equation, Equation (16), 

10 the first term is the constant k u which represents the voltage at the folly charged, zero current, steady 
state concentrations condition. It also represents the v 0 terms in the Nernst and Tafel solutions. The 
second term represents the nonlinear logarithmic dependence on current from the Tafel solution. The 
problem with the Tafel solution is that it assumes the current is large in one direction, but not both, 
while the approximation must be applicable to currents in both directions. Taking this into 

1 5 consideration yields the compromise term with the abs and sgn functions. The third and fourth terms 
represent the Nernst dependence on concentrations on either side of the charge storage reaction (and 
ignoring the fact that one reactant is in the solid state). The third term represents all the collective 
diffusion terms, as well as the electrolyte ions that are being consumed or generated at the electrode 
surface. As this concentration goes to zero, the terminal voltage gets small. As this concentration 

20 goes up, it has less and less effect on the voltage. The fourth and last term represents the effect of 
stored charge. As the stored charge approaches one, the voltage gets big. If it ever gets to zero, the 
voltage will also get small. 

Other major benefits to the linear-in-the-parameters approximate solution is that it is possible 
to identify the 4 ^-parameters directly from charge-discharge data using the batch or recursive least 

25 squares method. The requirement to do this is the actual measured terminal voltage and current, 
along with the estimated transient values of the stored charge and the diffusive material. 
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Note also that the charge transfer resistance through the double layer is contained 

implicitly in the current term in the linear-in-the-parameters approximation. Likewise, the double- 
layer capacitance explicitly is not modeled, as it is also implicit in the linear-in-the-parameters 
approximation as well as the diffusive term. 

5 

Charge Storage 

The charge storage will be through charge accumulation by integrating the current. Using 
the electrical analogy, this element is effectively capacitive. Any charge that passes into or out of the 
cell, will directly effect the concentration of the stored charge via Faraday's law, 
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c(0- Kt) . (19) 



Diffusion Behavior 

The diffusive components are lumped into one fractional-order term. This term effectively 
1 5 represents the diffusion of material toward or away from the electrode through the electrolyte, and 
with our assumed reaction, is thus on the same side of the reaction equation as the stored charge 



S + D + e reaction products. 

* charge 



20 As a fractional-order element, its transfer function is 



= k d *~* (20) 



where s is the Laplace variable, and and q are constant to be determined from data. However, the 
25 difiiisive material concentration must return to the assumed bulk value when the system is quiescent. 
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Thus the equation is converted to the time-domain and modified to reflect this return to bulk 



concentration 



cAO-cw-k,, 0 d;<m. (2i) 

5 

It should be remembered that fractional-order terms have time-varying initialization functions 
associated with them. These were determined to have an overall effect on the diffusive terms of less 
than 3% our specific study, and were not included (ignored). If the implemented battery model is 
used in other ways after this study, the initialization functions may need to be included, which might 
10 also require further adjustment of the constants in Equation (21). 



Self-discharge 

Self-discharge is implemented through a simple resistive (to flow) loss across the charge 
storage element, 

is - (22) 

K sd 

This should be thought of as a pathway for stored material to leak out of storage (bigger R s <t 9 less 
loss). 

20 Ohmic Cell Resistance 

Most conducting media provide some resistance to the passage of electrical current. One 
could attempt to use one of the many ways that this is often calculated in electrochemical cells for our 
model, but in this case the current term in the linear-in-the-parameters electrode equation is used to 
implicitly represent this resistance. If a value of resistance is needed, it can be obtained by 
25 linearizing the linear-in-the-parameters electrode equation. This is done here since the resistance of 
a cell is often discussed. Noting that resistance can be defined from Ohm's law as 
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(23) 



then 




(24) 



5 



or 



R. 



'ohmlc 



dv rf(Ar 2 ln(l^li|)sgn(/)) 
di di 



(24) 



One can evaluate tins for any steady-state current. Evaluating at zero steady-state current for 
symmetry gives 



Thus the linear-in-the-parameters electrode equation gives a direct measure of the total Ohmic 
resistance in the cell at zero current, which includes both the charge transfer resistance at the 
1 5 electrodes, and the Ohmic resistance of the electrolyte. 

Essen tialized Model Overview 

The model structure is given in Fig. 4. The system dynamic equations are repeated below to 
highlight the required parameters. 



10 




(25) 



20 



Terminal behavior: 



current into the battery is +/, 



(26) 



the terminal voltage is +v, 



stored material with self-discharge: 



(27) 
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diffusing material: 



c d (f)=c bulk -k d o^i'W, 



(28) 



electrode equation: v(c, , c d , 0 = k x + k 2 ln(l + 1 i |)sgn(0 + k 3 ]n(c d ) + k 4 In(l - c s ) .(29) 

5 

From this we can see that the required parameters are, self-discharge resistance R sd , fractional-order 
q and constant k di bulk concentration of the diffusing material c bulh , and the linear-in-the-parameters 
electrode equation constants, k u k 29 £ 3 , and 

1 0 Parameter Determination 

To most easily discuss the model parameter determination approach, it will be applied to a 
specific cell. However, it is believed that this model is equally applicable to other types of batteries 
and/or cells. In any event, the cell chosen is the NSWC Crane Pack ID 3602G. It is rated at 65 AHr, 
and uses 31% KOH concentration, and is maintained at 10 degrees C. The charge-discharge profile 

15 is a square wave, with 35% depth-of-discharge (DOD) and 104% recharge ratio. The square wave 
current is 26.29 A for 54 minutes charging, and -37.92 A for 36 minutes during discharge. Note that 
65 AHr= 3900 AMin, 35% of 65 AHr = 1365.1 AMin = 26.29A*54Min, 3900AMin- 1365.1 AMin 
= 2534.9 AMin, and 104% of 1365.1 AMin =1419.7 AMin = 37.92A*36Min. 

The self-discharge resistance is determined by the recharge ratio of 104%. It is assumed that 

20 roughly half of the extra 4% goes into self-discharge, and that the other 2% goes into overcharge 
reactions. The division of the excess could be determined more accurately using the Buder method, 
but that has not been done at this time. Assuming no current in or out, it is assumed that 2% of the 
stored material (1.3, AHr = 78 AMin) is lost during each 90 Min cycle via self discharge. Using this 
in the known solution to Equation (27) 

25 

c,(0=^(0)e" // ^ (30) 

gives 
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c s (/)=3900- 78 = 3900<T 90/ ** / . 
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(3D 



Solving this for the self-discharge resistance gives R sd =4455 ohms. This value was later modified 
after the simulation was implemented to force the constant-plus-taper charging to give a final charge 
5 of 3900AMin. The modified value of the self-discharge resistance that is used in the charging studies 
is R sd =4024 ohms, or equivalently 1/*^ =0.0002485 . 

To determine q and for the diffusive material, it is observed that most NiH 2 cells are 
diffusion limited. This must be true, since the terminal voltages get small after only 35% of the 
active stored material is used. To start, it is assumed that the bulk concentration is unity, 

10 



Thus we will assume that the diffusive material goes from unity concentration to roughly zero 
concentration as 35% of the stored material is used. With a constant current of -37.92 A, the solution 
15 to the diffusing material equation using Laplace transforms is 



s s g s 



The time response obtained by inverse transforming this is 

20 

c d {f)= 1 -37.92V*- < 34 ) 



To determine the unknown constants, we will evaluate this at the end of discharge, or £=36 Min, at 
which point the concentration goes to roughly zero, 

25 

l-37.92* rf (36)«=0, 

or 
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37.92* rf (36) 9 =l. 
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(35) 



Taking the natural logarithm of both sides and rearranging gives 

5 - ln(37.92) = q ln(36) + \n(k d ) . (36) 

This gives one equation for our two unknowns. If we now assume that the charging occurs similarly, 
we can obtain another equation. This is not exactly true, since the diffusive material will recover on 
its own, with no driving current, but this recovery is somewhat slower than the forced recovery. We 
10 will thus ignore the associated initialization response during the charging phase. Thus, charging at 
26.29 A for 54 Min yields the equation 

- ln(26.29) = q ln(54) + \n(Jc d ) . (37) 

15 These two equations can now be solved for the two unknowns, which become gr=0.9034 and 
Jfc<r=0.001036. 

To implement the diffusion term, an approximation to the fractional operator is necessary. 
The resulting approximation is 



20 



= 0.001036 0 01^^1 1444, + 0.3455 
s 2 +0.5078^ + 0.0013 



which is designed with a maximum Bode magnitude response error of 2dB between co — 10* 3 and 
co — 10 +2 . Better approximations containing more terms, or simpler approximations containing fewer 
25 terms are easily generated using the method Again it should be noted that the initialization terms 
have been ignored. 
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The parameters in the electrode equation are obtained by using a batch least squares 

optimization on the charge-discharge data from the actual cell. The resulting equation is 



v=. .» 



+. .« > ( +|j|)sgn(i)+0.02291n( C</ )~0.02621n(l.005*3900-c 5 ).(39) 



5 



Notice that the last term has been rescaled so that c a can go to the AHr capacity of the battery, which 
is 3900AMin. 

Model Verification Studies 

10 Figs. 5, 6, and 7 show actual charge-discharge data (the e o* points) from the cell discussed 

above. These figures also contain the charge-discharge plot (solid line) from the identified model 
also discussed above. The three plots are chosen at 500 cycles, 10000 cycles, and 40000 cycles 
respectively to shown how the specific cell has changed with time, and to show the ability of the 
model to track changes in the cell. 
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Damage Modeling in NiH 2 Batteries 



Background 

If the real-time observer is going to be able to provide information concerning the state of 
cell health, some understanding of the damage process in an electrochemical cell is necessary. There 

20 is no known instantaneous damage rate model for electrochemical systems, however, fatigue and 
creep damage modeling has been used for damage mitigating control in aerospace propulsion 
systems. In that application, damage per cycle information has been used, as well as a continuum 
damage model developed. The continuum damage model predicts the current damage rate as a 
function of local stress and strain, and can be used to implement a controller that trades a small 

25 amount of performance (dynamic response) for large gains in system lifetime. A continuum damage 
model will be required for work with batteries and fuel cells. 
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For electrochemical cells and batteries, a number of damage mechanisms have been 

identified which include: corrosion, electrolyte communication between cells, separator dry out, 

popping damage, plaque expansion and electrochemical degradation. Still other mechanisms may be 

present which have not yet been identified. However it still remains to quantify these mechanisms, 

5 determine their importance to a particular design, and to determine the operational relationships to 

damage accumulation. As these effects are quantified, the resulting relationships can be used to 

refine battery damage prediction. The next two sections discuss in a general way the mechanisms 

that may be responsible for battery cyclic damage. 

1 0 Damage Mechanisms for NiH 2 Batteries 

There are many ways to damage any cell. The most significant damage mechanisms for 
NiH 2 cells will be evaluated. 

A. Formation of y - phase NiOOH: 

15 At the positive nickel electrode, the charged material is NiOOH, and the discharged material 

is Ni(OH) 2 . As in many solids, however, these compounds can exists in difference solid phases. The 
normal charge-discharge cycling occurs between the beta phases, 

discharge ^ 

P -NiOOH ^-Ni(OH) 2 . 

^ charge 

20 

On overcharge, however, the excess heating due to 0 2 formation along with any large currents, 
causes the fi -NiOOH to undergo a phase change into the y -phase, y -NiOOH possesses two 
properties of interest to our studies. One is that it generates a lower voltage with respect to the KOH 
electrolyte, roughly IV and opposed to 1.5V for /? -NiOOH. The other physical property is that the 
25 spaces between the molecules are much larger for y -NiOOH than for J3 -NiOOH. Thus the 
formation of y -NiOOH can physically damage a cell stack due to the changing size associated with 
the phase change. Additionally, as the amount of y -NiOOH builds up, the usable charge storage. 
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capacity of the cell goes down, as the x -NiOOH cannot provide the voltage necessary for the 

electrical loads. 

Cells with a build-up of y -NiOOH can be reconditioned using the reaction mechanism 
proposed by Bode. This is shown in Fig. 8. It should be observed that the y -NiOOH will cycle 
5 electrically with a -Ni(OH) 2 . The original J3 -Ni(OH) 2 can be recovered by allowing the cell to 
stand for an extended period of time in a KOH solution. 

B. Formation of 0 2 : 

Overcharging a cell can occur by either continuing to charge the cell after all the -Ni(OH) 2 
10 has been converted to /? -NiOOH, or when the charging current is so large that the >0-Ni(OH) 2 
cannot react fast enough to form /3 -NiOOH. In either case, the excess current goes into the 
electrolysis of the cell electrolyte. The effect of this is the formation of 0 2 at the nickel electrode, 
along with heating. The problem is further complicated by the formation of oxygen in a hydrogen 
rich environment These can recombine chemically and release significant heat. This recombination 
15 can also do physical damage to the electrode. This chemical recombination of oxygen and hydrogen 
is often called popping. Clearly, overcharging is to be avoided. 

C. Heating: 

Heating is another major problem in NiH 2 batteries. As discussed above, it results in the 
20 formation of y -NiOOH, which reduces cell capacity and does physical damage to the cell. Another 
problem with heating is the increase in self-discharge reaction rates, which can also be manifested as 
a reduction in cell capacity. Any excess heat production must be avoided. 

Thermal Analysis and Hysteresis 

25 

This section will discuss the gross energy balances taking place in an electrochemical 
reaction. One needs to consider the free energy of the reaction, 
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AG = AH - TAS . (40) 



Here AG is the Gibbs free energy of the reaction and represents the energy available for external 
5 work, and in electrochemical cells is effectively the source of the terminal voltage. AH is the 
enthalpy of the reaction and represents to total chemical energy change during the reaction. The 
TAS term is the temperature multiplied by the entropy change of the reaction. The entropy change 
involves the reordering of the reacting molecules as they undergo chemical reaction. Multiplying the 
entropy change by the temperature of the reactants, TAS , gives the heat released or absorbed during 
10 the reaction. If the free energy is negative, the reaction will be spontaneous. 

Dividing the above equation by -wF, where n is the number of electrons involved in the 
reaction and F is Faraday's constant, gives 



AG _ AH } TAS (41) 
nF nF nF 

15 

The sign convention of a spontaneous reaction giving a positive voltage in an electrochemical cell 
then gives 

AG = -/zFv. (42) 

20 

Using this in the energy equation yields 



TAS 
nF 



25 Here 



AH = - nFv 



TN 



(43) 



(44) 
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is the thennoneutral voltage, or the voltage of the electrode reaction when there is no heat being 
generated Although free energies, entropies, and enthalpies are all generally dependent on the 
concentrations of reactants, the enthalpy of reaction for many electrochemical reactions is roughly 
5 constant. This allows the thermoneutral voltage to be used as a measure of when there is heating at 
the reaction site by comparing it with the actual terminal voltage. 

The energy flow in the reaction can be obtained by multiplying the above equation by the 

current, 

(45) 
(46) 



10 



. TAS 
nF 



or 



iv = iv m + Q r . 



Here Q r is the rate of heat release or absorption at the reaction site, and represents the unavailable 
1 5 electrical power. This equation can be rearranged to give 

a=»-(v-v m ) (47) 



which shows that comparing the terminal voltage with the thermoneutral voltage allows 
20 determination of the heating due to the electrochemical reaction. 



The above analysis has only considered the energy related to the electrochemical reaction, 
and has not considered any of the electrical behavior elsewhere in the overall cell, particularly any 
electrical losses. Whenever current flows in a conductor, there are electrical losses due to the heating 
25 of the material. In this case, the total Ohmic resistance as the primary electrical heating source will 
be used. As such, Equation (43) can be written as 
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v = v ro+ ^+//W. (48) 



This shows that the terminal voltage varies not only with the free energy of the reaction, but also with 
the resistive loss in the cell. Multiplying this equation by the current gives 

5 

TAS 

iv = iv m +i— + i 2 R 0hmfc . (49) 



This equation can be expressed in words as: 

10 electrical power = electrical power of the reaction + reaction heating/cooling + electrical 

heating. 

Thus the overall heating is a combination of electrical and chemical effects. 

15 The energy flow can be further generalized, by including the energy that goes into oxygen 

production during overcharge. Equation (49) is generalized adding this electrolyzing current, 

iv = Q-i 02 >m +0'-*o2>-^r + itnir-tXotm*) + '^oyc - ( 50 ) 

20 In words, this says: 

electrical power in = electrical power into the reaction + reaction heating/cooling 

+ heating due to O2 production + electrical heating. 

25 The discussion of this section can be summarized in Figs. 9 and 10 which show charge- 

discharge curves on a v-q diagram, along with the associated energy distribution as described by 
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Equation (49). These diagrams are somewhat ideal in the sense that they do not have self-discharge 

or overcharge effects included. In the case of overcharge, the figure would need to be modified to 

reflect the current going to the oxygen production as described in Equation (50). Another important 

effect that is not considered in these figures, or the discussion above, is that the voltage can move up 

5 or down as the diffusing material relaxes (independent of the stored charge). Since only the stored 

charge is plotted on the horizontal axis, the figures do not take the concentration of diffusing material 

into effect. 

At this point, sufficient data does not exist to quantify these effects for a particular battery 
design. Thus, based on available cycle life vs depth of discharge data an empirical continuum 
10 damage model for cells based on external variables, such as voltage and current, and internal variable 
such as species concentrations has been developed. It may be possible to explicitly include other 
internal variables, such as pressure and temperature in the future. This damage model allows the 
development of an instantaneous damage rate sensor and a damage mitigating control for NiH 2 cells 
which is discussed in the next section. 

15 

Battery Continuum Damage Modeling 

This section develops methods for using battery cyclic life data to infer instantaneous 
damage rate within the battery. This instantaneous damage rate is required to design optimum 
charging profiles and to establish control algorithms. Cycle life as a function of depth of discharge 
20 (DOD) data for NiH 2 cells is given for a temperature of 10 degrees C, as follows 



25 





Cycles 


DOD 


To 




Failure 


35% 


38,000 


50% 


19,000 
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5 For this data, a charge/discharge cycle has a 90-minute period, with a 54 minute charging 

duration. The recharge ratio is always less than 110%, and the recharge current was constant until 
94% recharge at which point a taper charge was used to finish the charging. The cycle life data is 
reproduced in Fig. 11. In this figure, only two data points are available at 10 degrees C. Four curves 
have been fit to the data, the Green-Hoffman model, a modified Green-Hoffman model discussed 
10 later, a power law, and a linear curve. Although only the two Green-Hoffman approaches are 
considered, any curve can be fit through two points. 

The Green-Hoffman approach uses the following exponential model 

N f GH =1885.04e 4 - 621(1 - DOD) , at 10° C (51) 

15 

where N f is the number of cycles to failure. It is assumed that when the total accumulated damage 
becomes equal to one, the cell has foiled. Otherwise stated, the total accumulated damage is 
normalized to one. Further assuming that the damage per cycle is constant for a given DOD, or more 
specifically, that the damage per cycle does not increase as the accumulated damage increases, then 
20 the damage per cycle is 

= 0.0005305*- 4 - 621 (1 - TCD > =— — . (52) 

N fCH 

Based on similar previous studies of mechanical fatigue, an expression for die instantaneous 
25 damage, D(f) , is derived based on the particular cycle life model. This is possible because the cycle 
life data implicitly contains all the information about how heating, overcharge, high current, and 
other hard to model factors damage a cell. The cyclic damage model expressed by Equations (51) 
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and (52) characterizes the damage hysteresis cycle size by DOD, which is actually representative of 

the charge amplitude of the cycle. The hysteresis cycle size could just as well be characterized by the 

voltage amplitude of the cycle. The purpose of the models to be derived here is to optimize the 

charging current profile. Models based on DOD as characterizing the size of the damage cycle are 

5 thus not useful because the overall charge to be stored during charging is always the same, no matter 

what charging profile is used. An alternative and more appropriate descriptor of cycle size for 

damage modeling is the voltage amplitude of the cycle. The voltage amplitude referenced to the 

minimum voltage point of the cycle (i.e. the discharge cutoff voltage (35%)) will be used to describe 

the damage cycle. The stress variable then will be the voltage difference based on this reference 

10 voltage. 

For a full charge/discharge cycle, the voltage amplitude v a varies with DOD (charge 
amplitude). This relationship is certainly a nonlinear function, but we will use a linear 
approximation. We express this as 

15 DOD=c x v a . (53) 

To obtain c x , it is adjusted until the simulated damage per cycle for a constant-plus-taper charge is 
equal to that given by Equation (52) for 35% DOD. This equation is substituted back into the 
damage per cycle, Equation (52), giving 

20 

D cyc = 0.0005305 e -*** M ™^° . (54) 

Now the accumulated damage as a charge/discharge cycle is traversed is the damage per cycle. Thus 
the integration of the instantaneous damage through the cycle is 

25 

\8(y)dv = = 0.0005305^ 62l+4621ClV * . (55) 

cycle 
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J(v) = ^. (56) 
dv 



5 It is now assumed that all of the damage is accumulated over the charging part of the cycle, that is the 
"uphill" part of the process, then 

p(v)rfv = = 0.0005305e- 462l+4621ClV - = 5.222 x KrV' 621 ^ . (57) 

v min 

10 It is noted that models including damage on the discharging portion of the cycle are readily possible 
using techniques similar to those used here. Changing the limits of integration gives 



v max — v min 

pO' + v^rfv = 0.0005305e^ 1+4621ciV - - 5.222 xlO' 6 e*' 62lc * v ° . (58) 
o 

15 Since v a = v max - , One can now infer an approximation to the shifted voltage referred damage 
rate, ^(v + v^), as 



S(y + v min ) = 0.0005305*- 4 621 4.621 cp 4 " 1 * v (59) 

20 The reason that this is approximate is due to the fact that at cyclic amplitude of DOD = v a = 0 , the 
Green-Hoffman model does not yield zero cyclic-damage, as we require. In fact, as discussed above, 
zero amplitude would actually allow healing. Then Equation (59) becomes 

S(y) =2.4129x 10- 5 c ie 462,Cl(v " Vrain) . (60) 
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This result in Equation (52) gives 

Va j3(y)dv = 5.222xl0^(e 462,civ " -l)= 5.222x 10" 6 e 462lCl v ° . (61) 

5 

Once the voltage referred damage rate is obtained, the instantaneous damage rate can be obtained 
by simple application of the chain rule 

w4=Tr^- (62) 

dt dv dt 

10 

The model requires that damage, D(t), be a monotonically increasing function, that is, the model will 
be conservative in that, healing of the battery will not be not allowed. Therefore the instantaneous 
damage rate is required to always be positive, thus Equation (62) is modified to 



15 *».§ .£ 

dt dv 



dv 



dt 



= <?(v(0)|v(0|. (63) 



Substituting the voltage referred damage rate into this equation gives the instantaneous damage rate 
as 

20 D(t) = 5.222 x 1 0^ c, e 4<a,Cl (v_Vnli " ) |v(/)| . (64) 

Choosing v^^X.lW, gives at last 

£>(f) = 5.222x 10-^c, e 4jaic>iv - 12) |v(0| • (65) 

25 
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This is the continuum damage model based on the Green-Hoffman cyclic damage model. 



Modified Damage Model 

The Green-Hoffinan exponential model, Equations (51) and (52), implies that cyclic damage 
5 occurs when the depth of discharge is zero. This is not physically plausible for an operational battery 
damage model. The degradation that occurs during shelf storage is not the same as cyclic damage. 
Thus, the Green-Hoffinan model is modified so that cyclic damage goes to zero as depth of discharge 
goes to zero. For the model to allow zero damage at zero DOD it is required that the cyclic damage 
have the form 

10 

D*. =c 2 (e^ D -e*°) = c 2 (e* MD (66) 

The constants c 2 and c 3 are determined numerically to best fit the data presented in the table at the 
beginning of this section. This process yields the damage per cycle as 

15 

= 1.0404 x 10- 5 ( e >*° 2D0D - 1) = J- . (67) 

For DOD =35%, then this equation gives £> ov = 2.6316xl0" 5 . To get the instantaneous damage, a 
linear relationship between DOD (the charge amplitude) and cyclic amplitude of voltage deviation is 
20 assumed, 

DOD=c x v a . (68) 

To obtain c l% it is adjusted until the simulated damage per cycle for a constant-plus-taper charge is 
25 equal to that given by Equation (67) for 35% DOD. Finding c, must wait until the end of the process 
however. This equation can be substituted back into the damage per cycle 



32 



WO 03/103071 



PCT/US03/17028 



= 1.0404x 10- s (e 3602 ™ - 1) (69) 

Then, the integral of voltage referred damage rate, <5(v) , for one charging cycle must be equal to the 
5 damage per cycle, and assuming damage formation only during the during the charging phase of the 
cycle 

V "p(v) dv = = 1 .0404 x 10" s (e 3 602ci v " - 1). (70) 

v niiD 

10 where 



S(v) = ^ (71) 
dv 



is the voltage referred damage rate and D (f) is the accumulated damage. Changing the limits of 
1 5 integration of Equation (70) gives 



v max "nun w a t \ 

J S(y + v min )dv = jSiv + v^dv = 1.0404X10" 5 (e 3 - 602 ™ -l) . (72) 

0 0 

The shifted voltage referred damage rate S(v + ) is now determined to be 

20 

<£(v + v min ) = l.0404xl0~ 5 (3.602c 1 )e 3 • 602c,, ' (73) 

or 

S(y) = 3.7475 x 10^ c, e 3 602ci ( - w » h) , v> VmiB . (74) 
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This voltage referred damage rate exactly satisfies Equation (70). Having the voltage referred 
damage rate, the instantaneous damage rate can again be obtained by application of the chain rule 



£>(/). f^ = ^ = J(v(0)v(f). (75) 
at dv dt 



Further requiring that the instantaneous damage rate is always positive gives 



dD dD 
D{t) = — 
dt dv 



dv 



dt 



- *(v(0) |v(0| . (76) 



10 

Substituting the incremental damage rate into this equation gives the instantaneous damage rate as 



D(t) m 3.7475 x 10" 4 c x e 3 ' 602 * {v - v ™> |v(/)| . (77) 

1 5 Choosing v min =l .2V, gives 

D(t) m 3.7475 x 10^ c x e 3602 *^ |v(/)| . (78) 

Now, adjusting c, until the simulated damage per cycle for a constant-plus-taper charge is equal to 
20 that given by Equation (67) for 35% DOD, that is D cyc = 2.6316xl0* 5 , yields c, =0.9707 . 

It should be noted that the instantaneous damage rate could also be expressed as a function of 
the current and rate of change of voltage with charge by continued application of the chain rule, 



. , N dD dD 



dv 
dq 



at\ 



dv 



dq 



41 



(79) 



25 
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Here dv/dq is the slope of the charging curve on the v-q diagram. 

In deriving Equation (76), it is required that the instantaneous damage rate always be 

positive. This was accomplished by using the absolute value of the voltage rate. This is justifiable in 

that whenever the voltage rate becomes large positive or negative, there are significant microscale 

5 changes taking place in the cell. An alternative result would be to have the damage rate go to zero 

whenever the voltage rate goes negative. This is justifiable in that the cell may not be undergoing 

damage when it is relaxing, that is, when the voltage rate is negative. Accordingly, the damage rate 

in this case becomes 



D(t) = 3.7475 x lO" 4 c x e 2 ^ {v ~ l 2) v{t\ v(/)2>0 (gQ) 
Z>(0 = 0, v(0<0 



In this case, the constant c, = 1.0901 . 

From the various forms for the damage, Equations (75)-(80), it can be seen that changing 
voltage is indicative of damage in the cell. The rate of change of voltage is an indication that the 
15 chemical reactions are proceeding in a manner to cause damage on a microscopic level. Thus 
keeping the voltage from changing abruptly should allow longer life in these cells. This would seem 
to indicate that a "soft take-off' and a "soft landing" are important in the charging process to avoid 
damaging the cell. Looking back at the energy balance equations with or without 0 2 production, 
Equations (49) or (50) give the following expression for the voltage rate 

nF at at 

In words this equation says: 

changing terminal voltage oc changing heat of reaction + changing heating 

due to current, 

25 that is, the heating rate is intrinsically tied to the damage process. Stated differently, or on a 
microscopic scale, the equation indicates 
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changing terminal voltage oc changing microscopic reordering rate + changing 



material flow rate. 

This is because the current is proportional to the concentration as indicated by Equations (9) 
5 and (10). Both of these processes involve the rate of change of material, that is, material 
motion and acceleration, which would be related to internal damage rate. This would seem 
to infer that damage rate and internal loss rates are related, thus justifying the use of voltage- 
difference rate as the continuum damage rate model stress variable. 

Now that an instantaneous damage rate model is available, one can use it to construct 
10 an instantaneous damage rate sensor and then to design an optimal charging profile that 
minimizes the accumulated damage per cycle. 

Effect of Damage in the Essentialized Model 

This section discusses the incorporation of the accumulated damage into the dynamic 
15 performance model. First it is important to associate the accumulated damage with some 
physical mechanism. The accumulated damage will be associated with the inability to store 
charge, or to accumulate the stored material. Thus this could be associated with the 
formation of the y - phase NiOOH, which actually does reduce the possible total stored 
charge. Likewise, overcharging always causes heating, which leads to the formation of y - 
20 phase NiOOH. Consequently, the accumulated damage will be associated with a reduction 
in the ability to accumulate the stored material. 

The inability to store charge could potentially be manifested into all of the parts of 
the performance model; Faraday's Law, the electrode equation (or the linear-in-the- 
parameters form), the charge storage term, the diffusion term, or the self-discharge. As the 
25 damage can often be associated with a reduction in the ability to accumulate the stored 
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material, a focus on incorporating the damage into the charge storage term will be 

considered. Equation (27) becomes 

stored material with self-discharge and damage: c s (/) = i(f) — c, (0 - ^(0 • (82) 

R sd 

5 But, as we have used the linear-in-the-parameters electrode equation, it must also be modified to 
indicate a loss of capacity. Thus Equation (29) becomes electrode equation with damage: 

v = k x +k 2 ln(l + | i |)sgn(i) + k 3 ln(c d ) + k 4 ln((l - £>(')) - c s ) (83) 
It will be appreciated that damage can be incorporated into the essentialized model in other ways, as 
is required by the specific device. 

10 

Damage Mitigating Control 

Control Philosophy 

Control system design is usually separated into at least two phases. Phase one is the 
15 design of optimal trajectories and associated inputs, that move a given plant from one 
operating condition to another, while minimizing some performance measure. This subject 
is the topic of the present section and the next Phase two is the design of a trajectory 
following controller (sometimes called a regulator or tracker) that provides a real-time 
control input perturbation to keep the plant operating near the designed optimal trajectory. 
20 This is discussed in the summary section. 

The preferred control philosophy is to charge the NiH2 cell in such a way that the 
damage incurred during the charging period is minimized, thus extending its cycle life. This 
is generally considered damage mitigating control or life-extending control. We will call 
die optimal charging profile life-extending charge. A damage mitigating control for use 
25 during discharge could also be considered, but the required discharge current was dependent 
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on time-varying mission needs, and would require hardware modifications that were 

considered outside the scope of this research. Some discussion of possible damage 

mitigating strategies for discharge is also included in the summary section. 



5 Performance Measure 

The performance measure to be minimized is the accumulated damage per recharge 
cycle. This can be obtained by integrating instantaneous damage rate with respect to time: 

J=\b(t)dt 9 (84) 

0 

where D(t) is the instantaneous damage rate which is obtained from either of the damage 
10 models or the damage rate sensor, Equations (65) or (78) , and t f for our problem is 54 Min. 
As it is, this performance measure can be minimized by using zero current, which does not 
recharge the cell. Thus, the performance measure must contain some information on the 
total charge required during charging. This is easily incorporated into a final state weighting 
term: 

15 J = Vfs(c s {t f )-c s (t f )J + /£>(/)<//, (85) 

0 

where c](t f ) is the desired stored charge at the end of charge ( <£l), w fs is the cost 

weighting associated with the error in stored charge at the end of charge, and again D(t) is 
the instantaneous damage rate which is obtained from either of the damage models or the 
damage rate sensor, Equations (65), (78), or (80), and t f for our problem is 54 Min. Note that 
20 only one weighting term is necessary as it can be either greater or less than unity, but is 
constrained to be positive. In the implementation, the self-discharge resistance is chosen to 
have a value which forces the final charge using constant-plus-taper charge (104% of 
1420AMin) to have a value of 3900AMin. The optimizer forces the life-extending charge to 
also have this final value. The constant charge (also 104% of 1420AMin) will have a final 
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charge larger than this since the self-discharge current is slightly less than the constant-plus- 



taper. 

Original Optimization Approach 

5 To determine the optimal charging profile, or life-extending charge, that minimizes 

the performance measure, Equation (85), a numerical optimization approach is used. The 
algorithm used is outlined in Fig. 12. Some typical charging results are shown in Figs. 13 A- 
D. 

1 0 Charging Profile Design for Specific Battery 

The optimization originally was accomplished by sampling the charging current in 
time, and optimizing each little sample increment. This optimization approach is very slow 
and requires many hours of computation time. The results shown in the last section, 
however, seemed to indicate that the optimal charging profile was similar to a parabola. 

15 Using this idea allows a much faster optimization to be accomplished, although the results 
are probably a little conservative as the truly optimal currents at each point in time are not 
being calculated. The fast optimization approach is outlined in Fig. 14. Here a ten point 
cubic spline is chosen to represent the current profile so as to allow some variation from a 
parabola (other time functions were also considered which gave similar results). Charging 

20 profiles are shown in the next several plots. In performing the optimization, it has been 
important to be sure to include the cost of requiring the zeroing of a non-zero current at the 
end of charging, as the large voltage rate there may incur significant damage. Also notice 
the constant current charge is predicted to last for only 25465 cycles, while the actual cell is 
still cycling after 50000 cycles. It is believed that this indicates either some problem with 

25 the cycle life data, or that this particular cell has used a design different from that of the 
cycle life data. 
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Recognizing that the cell voltage response is diffusion limited on discharge, along 

with the knowledge that much of the damage occurs near the end of charge when most of the 

overcharge occurs, allows one to try a charge-discharge profile that backs off from the 100% 

charged point Said otherwise, the life-extending charge profile from 60% to 95% fully 

5 charged can be determined. On discharge, the cell will still use 35% of the AHr capacity 

since it is diffusion limited at cutoff. On charge, the cell will be kept away from the 

overcharging region, and thus reduce the damage per cycle. The resulting life-extending 

control will be implementable in real-time due to the use of an on-board dynamic observer 

which can keep track of the stored charge material. The use of a real-time observer is 

1 0 discussed in the next section. 

The expected improvement associated with the life-extending charge is summarized 
in Tables I, II, and m which show a significant potential extension in life expectancy. The 
percentage life increase should be regarded as approximate due to the use of the fast 
optimizer and relatively large simulation timestep. It is interesting that the optimal solution 

15 using the \dvldt\ term in the damage rate model (Table I and Figs. 15 A-F), gives the same 
optimal charging curve as incurring damage only when dv/dt is positive (Table III and Figs. 
17 A-F), although the actual cycle life is different because of different damage rate models. 



Damage 


Cycles to 


Life-extending-Charge 


%Life 


per Cycle 


Failure 






Extension 






Damage per 


Cycles to 








65-100% Cycle 


Failure 
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Constant- 
Current 

Charge 


3.9269e-005 


25465 


2.0780e-005 


48123 


88.98% 


Constant + 
Taper 

Charge 


2.6316e-005 


38000 


2.0780e-005 


48123 


26.64% 



Table I Life-extending charge compared to constant current and constant current + taper. 





Damage 
per Cycle 


Cycles to 
Failure 


Life-extending-Charge 


% Life 
Extension 


Damage per 
60-95% Cycle 


Cycles to 
Failure 


Constant- 
Current 
Charge 


3.9269e-005 


25465 


1.9535e-005 


51190 


101.02% 


Constant + 
Taper 

Charge 


2.6316e-005 


38000 


1.9535e-005 


51190 


34.71% 



5 Table II 60%-95% Life-extending charge compared to constant current and constant current + 

taper. 
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Damage 
per Cycle 


Cycles to 
Failure 


Life-extending-Charge 


%Life 
Extension 


Damage per 
65-100% Cycle 


Cycles to 
Failure 


Constant- 
Current 
Charge 


3.9269e-005 


25465 


2.2080e-005 


45290 


77.85% 


Constant + 
Taper 

Charge 


2.6316e-005 


38000 


2.2080e-005 


45290 


19.18% 



Table HE Life-extending charge compared to constant current and constant current + taper 

using only positive dv/dt in the damage rate law.. 



5 Suggested Real-time Implementation 



Real-time Observers 

An important result of state space control theory is the concept of a state observer. 
The total number of dynamic "states," or "state variables," in a system is equal to the number 
10 of energy storage elements in the system, which is also equivalent to the number of time 
derivatives that are needed in the differential equations representing that system. The "state" 
of a given system is the collection of all the "state variables" of that system at a given point 
in time. This collection is all the information needed to advance the "state variables" 
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forward in time through the differential equations representing the system. With these 

definitions, the concept of a state observer will now be clarified. 

A state observer is effectively a dynamic model of a given system that runs in real- 
time, and in parallel with the actual system. Typically, only a small number of the dynamic 
5 states of an actual system are measurable, or are sensed. In the case of an electrochemical 
cell, typically only terminal voltages and currents, and perhaps some pressures and 
temperatures are all that are sensed. The state observer is then used to provide an estimate of 
the non-sensed variables. There are usually two forms for state observers, the open-loop 
observer and the closed-loop observer. 

10 An open-loop observer receives the same inputs as fee actual system, and then 

provides the collection of dynamic states as an output. If the model being used in the 
observer is accurate, then all the dynamic states from the observer will accurately represent 
the true states in the actual system. This has been applied successfully to jet engines, as well 
as other places. A closed-loop observer receives the same inputs as the actual system, as 

15 well as all of the sensed dynamic states and outputs from the actual system. Comparing 
these sensed variables with the corresponding simulated variables, the difference can be used 
as a feedback signal to force the observer dynamic states to converge to the true dynamic 
states of the actual system (see Fig. 2). This approach does not require the model to be as 
accurate as in the open-loop observer due the correction of the feedback from the sensed 

20 variables. Closed-loop observers are used in most state space control designs, where they 
have also been used to filter noise out of the measurements. 

An optimal controller for an electrochemical cell as shown in Fig. 18 is designated 
generally by the numeral 100 and incorporates the hybrid cell model, as discussed 
previously, into a microprocessor based simulation that is run in real-time. Sensed variables 

25 from an actual cell 102, which could be any rechargeable battery device, are input to a 
simulation microprocessor 104. The internal states of an exemplary cell, include, but are not 
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limited to normal charge-discharge operation, and overcharge reactions. These internal 

states were previously discussed and are embodied in Equations 1-6. The processor 104 

receives inputs such as voltage, current, temperature, and any other quantifiable state 

variable in the form of an input signal 106, as well as known system parameters such as 

5 charge storage properties of the cell, the diffusive behavior of the cell, charge transfer 

resistance, capacitance, self-discharge resistance, and the like. Likewise, a measured output 

signal 108 representing the comparable state variables mentioned above are generated by the 

cell 102. Both the input signal 106 and the output signal 108 are received by the 

microprocessor 104. As appreciated by the skilled artisan, the microprocessor incorporated 

10 into the microprocessor 104 contains the necessary hardware, software and memory for 
implementing the teachings of the present invention. 

The microprocessor 104 includes a feedback component 1 10 and a real-time observer 
component 112. Generally, the sensor 104 forces the dynamic states of the real-time cell 
observer 112 (simulation) to converge to the corresponding states in the actual cell 102. The 

15 feedback component 1 10 is crucial for this problem as the actual system is highly nonlinear 
as well as spatially distributed, and the model dynamic states represent aggregations of many 
actual states. Even so, the feedback to the observer will force the errors in the measured 
variables to zero, and the unmeasured variables will accurately represent the corresponding 
actual variables. The equations associated with the feedback process can be set out as 

20 New states estimates = old state estimates + Gain x Error signal. 

And a similar equation holds for the parameter identification 

New system parameters = Old state system parameters + Gain x Error signal. 
These equations are updated every sample time in the microprocessor 104. The equations 
have also been developed for fractional-order system components. Equations 27 and 30-37 

25 are exemplary of the processes implemented by the feedback component 110. In particular, 
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the feedback component generates a correction signal 114 that represents a real-time estimate 

of the amount of damage being done to the cell during re-charging. 

The correction is received by the real-time observer component 112 which also 

receives the input signal 106. With these inputs an estimated output signal 1 16 is generated 

5 and returned to the feedback component 110 for converging the dynamic states of the 

observer component 1 12 to the corresponding states of the cell 102. As a result, the observer 

component 1 12 generates an estimated internal state signal 118. The estimated output signal 

116 provides a charging profile that minimizes damage during charging which, in turn, 

prolongs the life of the cell. It is expected that a real-time observer for an electrochemical 

10 cell will be able to provide accurate information about state-of-charge, active species 
concentrations, and internal temperature and pressure of the actual cell 102. The 
information included in the state signal 118 provides important information about the 
behavioral properties of the cell 102 under various charging and discharging conditions. The 
model structure employed by the observer component is shown in Fig. 4. Equations 26-29 

15 are exemplary of the process implemented by the observer component 112. Any theory 
required for the fractional-order components in the observer has already been developed. 
More advanced observer designs may also provide some indication of damage and changing 
system parameters as discussed below. 

20 Observer Based Control 

Typically real-time observers are used in two ways, for continuous-state feedback 
and for decision making. Much of modem control theory is based on state feedback. State 
feedback requires knowledge of all the dynamic states of the given system, which are then 
fed back to the system inputs to improve overall system performance. As most of these 

25 dynamic states are not measurable, most state feedback controllers use the estimated states 
from real-time observers for feedback purposes. 
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Observers can also be used as part of the general decision making process. To do so, 

then requires an extra layer of intelligence. For example, if the observer predicts that certain 

dynamic states have varied too far from where they are supposed to be, the higher level of 

intelligence can determine if some sensors have failed, if some actuators have failed, or if the 

5 actual system has become damaged in some way. Of course, this intelligent layer must be 
programmed by a knowledgeable user. Another decision making application could be to 
shut the system down if certain observed dynamic states exceeded those which would 
normally be deemed acceptable. These decision making control strategies will be addressed 
later when discussing advanced controllers and life-extending control. 

10 Some of the issues associated with observer based control of batteries and fuel cells 

require a clear understanding of improved performance. Some potential definitions of 
improved performance might include maximum energy out, maximum power out, constant 
voltage out, constant current out, constant power out, prolonging state of health, or minimum 
damage rate, among others. Other issues are more hardware related. For example, in a 

15 collection of cells, control of each cell would probably be important. An array of observers, 
one for each cell, would then be able to provide information about individual cells. As long 
as each cell could be independently controlled, this information could actually be used to 
handle any degradation or failures of individual cells. Of course this also requires the usage 
of power electronics to control the charging and discharging of individual cells. This could 

20 include the use of buck-boost DC-DC converters to provide constant voltages or currents 
out, as cell voltages and states of charge change. Finally, any theory necessary for control of 
the fractional order components will be developed. 

Advanced Real-time Observers 
25 An advanced optimal controller is designated by the numeral 200 and incorporates 

many if not all of the same teachings as the system 100. The controller 200 for 
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electrochemical cells, such as shown in Fig. 19, requires an extra level of intelligence over 

the standard real-time observer system. The system 200 is a microprocessor based system 

with additional components that facilitate optimal monitoring and re-charging of a galvanic 

cell. As in the system 100, the system 200 includes an actual cell (or multiple cells) 202. 

5 Sensed variables from the cell are input to a simulation microprocessor 204 which receives 

inputs such as voltage, current, temperature, and an other quantifiable state variables in the 

form of an input signal 206, as well as known system parameters such as the charge storage 

properties of the cell, the diffusive behavior of the cell, charge transfer resistance, 

capacitance, self-discharge resistance and fee like. The cell 202 generates a measured output 

10 signal 208 that is received by the microprocessor 204 which contains the necessary, 

hardware, software and memory for implementing the teachings of the present invention. 

The microprocessor 204 operates in much the same manner as the microprocessor 

104, and includes a feedback component 210 which generates a correction signal 214 that is 

received by an advanced real-time observer 212. The correction signal 214 contains much 

15 the same type of information carried by signal 1 14 as discussed above. As in the previous 

embodiment 100, the microprocessor 204 receives the input signal 206 and the output signal 

208. The observer 212 generates an estimated output signal 216 which is returned to the 

feedback component 210, and an estimated internal state signal 218. The observer 212 is 

distinguishable from the observer 1 12 in that the cell model parameters contained therein are 

20 updated as they change with age in order to provide more accurate real-time estimates of the 

cell's internal states. In other words, by accurately modeling the internal parameters of the 

cell at different stages of the cell's lifetime, the observer 212 more closely tracks the 

performance of the cell to ensure that the re-charging cycle minimizes damage to the cell. 

In addition to the added features of the observer 212, the controller 200 includes a 

25 battery control system 220, a supervisory intelligent controller (SIC) 222 and an 

instantaneous damage rate sensor 223. These additional components facilitate the optimal 
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re-charging of the cell 202. The sensor 223 receives the input signal 206, the estimated 

output signal 216, and the estimated internal state signal 218. From these inputs,-the sensor 

223 generates a damage rate signal 224 which provides a real-time estimate of the amount of 

damage being done to the cell 202 during the re-charging process. The sensor 223 

5 implements the teachings embodied at least in Equations 61-65 discussed above. In any 

event, the damage rate signal 224 is sent to the SIC 222 and the battery charge control 220. 

In particular, the SIC 220 receives the damage rate 224, the estimated output signal 

216, the estimated internal state signal 218, and a desired performance input signal 226. The 

signal 226 defines such things as maximum desired power after charging, the time allowed 

10 for charging, desired life time performance of the cell and related performance data. These 
inputs 216, 218, 224 and 226 are processed by the SIC 222 to generate a charging profile 
signal 228 that is received by the battery charge control 220. The signal 228 effectively uses 
the desired performance of the cell to modify the input signal 206 generated by the control 
220 so as to trade-off performance of cell for the desired purpose of extending the life of the 

15 cell. In other words, the battery control system generates the input signal 206 in a manner 
that is continually updated to provide an optimal charging profile that does the least amount 
of damage over the life o the cell. The performance measure to be minimized is the 
accumulated damage per recharge cycle which is exemplified in Equation 85. Accordingly, 
the algorithm used for the identification of the cell hybrid-model parameters is preferably 

20 automated and is expressly set forth in Fig. 12. This allows two things. First, it allows the 
tracking of parameters as a cell changes during its lifetime. Some knowledge of parameter 
range could indicate when a cell was nearing failure. Second, it allows the on-line 
determination of the cell model parameters to begin with, or effectively, the correction of the 
initial parameters provided to the real-time observer 212. Effectively, the real-time observer 

25 212 could then 'learn" any cell. Preferably, the advanced controller 200 also includes 
models of any damage mechanisms as a function of the system states. And it is believed that 
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degradation algorithms should be included into an observer structure to give a real-time 

damage indicator for a given cell. 

The advanced controller 200 for an electrochemical cell would also require an extra 

level of intelligence as provided by the SIC 222. Specifically, the advanced controller would 

5 use the information provided by the advanced real-time observer 212 to maintain cell 
performance as the cell's parameters change with time. This would require a continuous 
change of controller gains as the effect of physical changes in a given cell were updated by 
the advanced real-time observer. 

Secondly, the advanced controller 200 would need to implement some of the decision 

10 making skills as discussed earlier in the Observer Based Control section. If the estimates 
provided by the real-time observer vary too much from the actual cell outputs, then decisions 
would be required concerning sensor, actuator, or cell health. If any failures were detected, 
additional intelligence would be required to accommodate a detected failure. Of course, the 
chosen battery and control structure must allow accommodation. 

15 Finally, assuming some form of a damage model was available through the use of a 

damage rate sensor, the controller could then make performance-lifetime trade-off decisions. 
Associated with this intelligence, there may need to be off-line optimization studies 
performed that determine the most efficient charging (or discharging) profile which will also 
prolong the cell lifetime. Ultimately, designing a cell charging (or discharging) profile that 

20 extends the useful life of the cell by minimizing accumulated damage would be a useful 
result. 

Summary and Further Considerations 

Summary 

A methodology has been presented for extending the life of NiH 2 cells via a damage- 
25 mitigating charging current This has required the development of an essentialized dynamic 
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model that will accurately reconstruct internal states in a given cell. The parameters of this 

model are readily obtainable from either off-line data, or adaptively on-line from real-time 

data. To determine the life-extending charging current, an instantaneous damage rate model 

that determines the damage rate at any point in time based on voltage measurements, has 

5 been developed. Using the instantaneous damage rate model with the essentialized cell 

model, allows the development of a damage mitigating charge profile for extending the life 

of the cell. This charge profile is obtained via optimization, and can be modified by trading 

recharge for damage. Once the life-extending charge profile has been obtained, the 

implementation of this profile is performed via a real-time controller that uses a real-time 

10 observer for reconstructing the internal cell variables based on current and voltage 

measurements at the cell terminals. This real-time controller implementation can be used 

over the entire life of the cell by updating any combination of the model parameters, control 

parameters, or charging profile based on the data from the present charge-discharge cycle. 



15 Model Development 

1. Developed an essentialized Nfflb cell model structure. 

2. Developed a method for determine the essentialized model parameters for charge- 
discharge data. 

3. Developed a deeper understanding of the charge-discharge hysteresis behavior in 
20 the voltage-charge plane. 

Damage Model 

1. Developed an instantaneous damage rate sensor for NiKfe cells based on cell 
voltage measurements. 

25 2. The instantaneous damage rate is shown to depend on voltage and voltage rate. 
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It is shown that the size of the voltage-charge hysteresis curve is related to the 



damage per cycle. 

Life-extending Charge Development 
5 1 . A life-extending charging profile has been developed by optimizing the charging 

current with respect to the damage per cycle and total recharge. 

2. This life-extending current should provide significant life extension over the 
present recharging methods. 

3. Significantly more life extension for a given cell can be obtained by not 
10 recharging the cell to 100%. 

4. Some heuristic rules for minimum damage charging that have emerged are: 

a. Increase the charging current slowly at the beginning of charge. 

b. Put more of the charge in during the first half of charging, 

c. Decrease the current to zero slowly at the end of charge. 

15 d. During charge, try not to allow the terminal voltage to increase much beyond 

the thermoneutral voltage. 

e. Try to keep the voltage-current hysteresis-curve area as small as possible. 



Observer Development 

20 1. A fractional-order observer theory (necessary for control implementation of the 

model) has been developed. 
2. A parameter identification technique is discussed for allowing a real-time 
essentialed model to track changes in its parameters in real-time over the life of a 
cell. 

25 3. A tracking controller for implementing fee optimal charging profile in real-time 

based on observed system states is discussed. 
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4. A learning controller that continuously updates the charging profile based on 

changing cell parameters has been discussed. 



Further Considerations 

5 

Overcharge Control 

It is noted that Eb pressure seems to be related to overcharge, and thus to damage. 
By monitoring the H2 pressure in a Nfflb cell, it is possible to provide some control for 
overcharging. The cell pressure was not included in the essentialized model, but could be 
10 without much increase in complexity. H2 pressure can then be related to both overcharge 
and damage. 

By using the instantaneous damage rate model, both a real-time damage rate sensor 
and an on-line damage monitor can be used. The damage rate sensor requires knowledge of 
the amount of stored material, and the present voltage. The damage monitor, or life 
15 predictor, requires the damage rate sensor, or can be implemented using an accelerated 
degradation model. 

Ultimately a charge-discharge controller can use feedback from these sensors to 
trade-off lifetime with performance. 



20 Other Chemistries and Fuel Cells 

It is expected that this methodology can be readily applied to other chemistries, such 
as Li-based cells if cycle-life data is available. If this data is not available, a cycle-life model 
can be inferred from an accelerated degradation approach. 

With regard to rechargeable fuel cells, this methodology should also be applicable, 
25 however it will be necessary to modify the essentialized model structure for the two fuel 
sources. It will also be necessary to understand the damage mechanisms associated with fuel 
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cells. Alternatively, it is expected that a tracking controller should be able to immediately 



improve performance in a fuel cell. 

Discharge Damage Control 

5 Although discharge control was considered outside the scope of this research, various 

possibilities exist to reduce damage during the discharge part of the cycle. It was assumed in 
the damage modeling that the entire damage during the cycle took place over the charging 
part of the cycle. It is of course possible to apportion that damage over both the charge in the 
discharge parts of the cycle. In the former, the damage model is dictated by the form of the 

10 damage per cycle and would remain very much the same. Future research has generalized 
the damage model to include damage on both charging and discharging. With the 
availability of the generalized damage model, various strategies can be investigated which 
would control damage during discharge. The additional aspect of such research is that it 
could influence the manner in which the battery power is allocated. The acceptability of this 

1 5 would depend on the specific application. 

Damage Testing and Future Data Needs 

The present research allowed access to the Crane database. From our observation of 
the database, which is certainly not a detailed analysis of the database, we have found it to be 

20 an excellent resource for evaluation of battery designs and battery chemistries, etc. 
However, it was not clear that appropriate data was available in the database, for evaluating 
and creating damage models for particular battery types. It has been our observation, that in 
some of the battery literature, the battery design aspects and the battery operational aspects 
are not clearly delineated. It is important to separate these two areas, and to develop as a 

25 separate discipline, battery operational theory. In particular, clean operational data seems to 
be required. By this we mean tests to failure, with key parameters held constant One 
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direction of future research would be to establish testing methods specific to obtaining usable 

damage (operational) models. This involves accurate assessment of the equation forms for 

cycle damage as well as proper statistical treatment of the data. 

Thus, it can be seen that the objects of the invention have been satisfied by the 

5 structure and its method for use presented above. While in accordance with the Patent 

Statutes, only the best mode and preferred embodiment has been presented and described in 

detail, it is to be understood that the invention is not limited thereto or thereby. Accordingly, 

for an appreciation of the true scope and breadth of the invention, reference should be made 

to the following claims. 
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